Method for Simulating the Movement of Particles 



BACKGROUND OF THE INVENTION 



[0001] The invention relates to a method for simulating the 
movement or distribution of particles in a medium under the 
influence of a changing phase interface between two 
neighboring phases in said medium and to a computer program 
product comprising software code to perform said simulation. 
A changing interface in the context of the invention means an 
interface which moves at least partially in the course of 
time . 



[0002] The properties of most semiconductor devices, such as 
MOSFETs , are determined by the spatial arrangement of so- 
called dopant or (impurity) atoms, such as e.g. boron or 
arsenic atoms. In the absence of dopant atoms, the conduc- 
tivity of the semiconductor is guite low. Regions with a high 
concentration of so called donor atoms (such as arsenic or 
phosphorus) are called n-type regions, because in such regions 
the concentration of negatively charged free electrons in the 
semiconductor is increased. Regions with a high concentration 
of so called acceptor atoms (such as boron or indium) are 
called p-type regions, because in such regions the concentra- 
tion of positively charged free holes is increased. A posi- 
tively charged hole can be understood as a missing electron. 



[0003] In the fabrication process of semiconductor devices, 
such as transistors, one of the main tasks is to achieve a 
certain spatial distribution of dopant atoms in the device. 
Usually, at the beginning of a fabrication process, one starts 
with a semiconductor wafer with a uniform dopant concentra- 
tion, and at the end of the fabrication process, one has 



1 



different doping concentrations in different regions of the 
device . 

[0004] In state-of-the art semiconductor technology, the 
following method is frequently used to fabricate regions with 
different concentration of dopant atoms: 

[0005] In a first step, dopant atoms are implanted into the 
semiconductor crystal by so-called 'ion implantation', i.e. by 
shooting dopant atoms at the semiconductor with high energy. 
Using different ion acceleration energies, one can adjust the 
average depth from the semiconductor surface, at which the 
ions come at rest in the semiconductor. By using masks one can 
achieve that the implanted ions will enter the semiconductor 
only at desired regions. 

[0006] In a second so-called annealing step, the semiconduc- 
tor wafer is subject to a high temperature treatment, e.g. 
treating the wafer 10 seconds at 1000°C, or 1 hour at 70°C. 
This is necessary to heal crystal damage which is caused by 
ion implantation, and to electrically activate the dopant 
atoms, i.e. to put them into semiconductor lattice sites, 
which requires some thermal energy. The temperatures needed to 
heal the crystal damage and to activate the dopant atoms are 
so high, that the dopant atoms diffuse in the semiconductor 
crystal. Therefore, not only the crystal damage and the 
activation of the dopant atoms changes, but also the spatial 
distribution of dopant atoms may change during an annealing 
step . 

[0007] Fabrication processes of semiconductor devices can be 
modeled with the help of computer simulations which are based 
on physical models, so called process simulations. A main goal 
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of process simulations is to calculate accurately the distri- 
bution and the electrical activation of dopant atoms in the 
completed semiconductor device. This allows to compute the 
electrical characteristics of the devices and to study the 
effect of changes in the fabrication process on the device 
characteristics. Such process simulations are done by most of 
the manufacturers of integrated circuits. An example for such 
a process simulation can be found in the European Patent 
Application EP 0 864 991. 

[0008] The implantation of a high dose of dopant atoms causes 
a lot of crystal damage to the exposed semiconductor surface. 
The reason is, that the incident high energy ions can collide 
with the semiconductor atoms, before they eventually come to 
rest, and kick semiconductor atoms out of their ideal lattice 
site. If the implantation dose is so high, that most of the 
semiconductor atoms in a certain volume of the silicon are 
kicked to some other position in the crystal, then the 
crystalline order of this volume is finally completely lost, 
and the region is considered to be amorphous. Those parts of 
the semiconductor wafer, which are not exposed to the ion 
bombardment or which are too far from the surface to be 
reached by the implanted ions, remain crystalline. 



[0009] Upon high temperature treatment, an amorphized region 
of semiconductor will recrystallize . If the amorphous region 
is neighboring a crystalline region, the crystallization 
process usually occurs in the form of a solid phase epitaxial 
(SPE) regrowth (see e.g. G. L. Olson and J. A. Roth, 'Kinetics 
of Solid Phase Crystallization in Amorphous Silicon', Mat. 
Sci. Rep. Vol 3, Nr. 1, June 1988, pp. 1-78). 

[0010] In SPE regrowth, the amorphous material crystallizes 
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along the interface between crystalline and amorphous mate- 
rial. The crystalline region acts as a seed for crystalliza- 
tion. Starting at the interface between the amorphous and the 
crystalline phase, the atoms of the amorphous part are 
rearranged in such a way, that they extend the crystalline 
microstructure of the crystalline part. By this mechanism, the 
crystalline region becomes larger and the amorphous region 
becomes smaller. The interface between amorphous and 
crystalline structure moves towards the amorphous part of the 
material, until eventually the complete amorphous part is 
recrystallized . 



[0011] The distribution of dopant atoms in the semiconductor 
can change in the course of time, e.g. during the SPE 
regrowth, either by diffusion, or by some force pushing the 
atoms in a certain direction. In particular, during SPE 
regrowth, dopant atoms, which are in the amorphous part of the 
device before the SPE, can be pushed in the direction of the 
recrystallization front. The driving force for this dopant- 
pushing is, that dopant atoms are in an energetically more 
favorable state in the amorphous region of the semiconductor. 
Near the crystallization front of the SPE, such dopant atoms 
prefer to stay at the amorphous side of the amorphous- 
crystalline interface, rather than being incorporated into the 
crystal structure of the crystalline part. In this situation, 
in average, each dopant atom which is originally situated in 
the amorphous part of the device may be shifted a certain 
distance until it is eventually built into the crystalline 
material . 



[0012] This dopant redistribution during SPE can have a 
significant impact on the final distribution and electrical 
activation of dopants in semiconductor devices. Therefore, it 
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has a significant impact on the electrical characteristics of 
such devices. 



[0013] Until now there is no process simulation available 
with which it is possible to determine accurately the movement 
and redistribution of the dopants during annealing steps in a 
medium having two different phases, particularly during SPE 
regrowth. As annealing is very essential to the production of 
semiconductor devices, it is of a great interest to determine 
the movement or distribution of the dopants by means of a 
computer simulation . 



SUMMARY OF THE INVENTION 



[0014] It is an object of the present invention to overcome 
this lack of process simulation work. This and other objects 
are solved with a method for the simulation of the particle 
movement and a computer program product for performing said 
method, as described below. 



[0015] According to the invention a method for simulating the 
particle movement at the interface between two phases of a 
medium is presented. This method is applicable to many 
problems involving movement of particles neighboring a phase 
interface, in particular it can be used to describe the dopant 
movement and distribution during and after SPE regrowth. 



[0016] The method is applicable on any specified particles, 
e.g. impurities, in any medium consisting of a material having 
at least two neighboring phases. These phases may be e.g. 
different crystallization patterns, but the method would also 
be applicable to e.g. solid and liquid phases. 
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[0017] The method comprises at least two steps for 
determining the movement and optionally at least a third step 
for determining the distribution of the particles. 

[0018] In a first step, the temporal and/or local evolution 
of the phase interface is determined. 

[0019] In a second step, the movement of the particles is 
determined in dependence of the previously determined 
evolution of the phase interface. 

[0020] In a preferred embodiment, the distribution of the 
particles at certain times during the movement is determined 
in a third step. 

[0021] Typically, these steps are automatically performed by 
means of a suitably programmed computer. 

[0022] Particularly, the distribution can be determined at 
the steady state position of the phase interface, which means 
that the phase interface remains constant in time and place, 
e.g. at the end of an annealing treatment. In the borderline 
case the steady state is reached only when the phase interface 
moves out of the medium which means that one of the phases 
disappears. The distribution of the particles can be 
preferably determined even in this borderline case. 

[0023] The phase interface can be described in two different 
models: Either there is a smooth transition between the two 
phases and the interface is extended over a non zero range in 
the medium, or the transition is discontinuous and the 
interface between the phases is sharp. According to the 
invention, it is possible to determine the movement or the 
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distribution of the particles on the basis of both models. 

[0024] The medium and in particular the evolution of the 
phases can be described by a data field wherein a first value 
represents the first phase and a second value the second 
phase, and values in between describe the transition. This 
data field can be representing a characteristic value of the 
medium, e.g. if one phase is crystalline that value can be the 
crystallinity or a related parameter. 

[0025] In a preferred embodiment, the determination of the 
temporal and/or local distribution of the phase interface can 
be calculated by an equation, in particular by a partial 
differential equation determining said data field. Assuming, 
that the first value for the first phase is 1 and the second 
value for the second phase is 0, a suitable differential 
equation is 

^ = div(agrad{C)) + pC\\-Cy(C-y) 
at 

wherein C is the characteristic value, div and grad are 
spatial differential operators and a, p and y are suitable 
parameters. The first term of this equation, div (a • grad ( C ) ) is 
a diffusion term. One can interpret this term in the sense, 
that the first phase diffuses (spreads out) a bit into the 
second phase and vice versa, with a being the diffusion 
coefficient. The second term of this equation, p*C* ( 1-C) * (C-y) , 
describes the trend, that in absence of any gradient of the 
characteristic value, a region, which is close to phase 1 
(l>C>y) turns completely into phase 1 (C converges to 1) , and 
a region which is close to phase 2 (0<C<y) turns completely 
into phase 2. The parameter P describes, how fast the final 



state (C=l or C=0) is reached, and the parameter y is the 
threshold value. The second term can be seen as an antipode of 
the first term, trying to establish a sharp interface between 
the two phases. 



[0026] For y= H, the combined effect of both terms in the 
equation would be, that the interface between the two phases 
is smeared out a bit and does not move. For y < ^, the 
interface would move in direction of phase 2, because the 
integral over the termp^C* (1-C) * (C-y) has a positive value in 
the transition region. 

[0027] Generally, the parameters a, P, y determine the width of 
the transition region between the phases and the time scale of 
the phase change at the interface. 

[0028] The data field can have an influence on other data fields, 
e.g. the particle distribution can be described by a data field, the 
time evolution of which is influenced by the data field describing 
the characteristic value for the phases. 

[0029] It is possible to determine the evolution of the phase 
interface by the calculation of probabilities for regions adjacent 
to the phase interface to switch from one phase into the other in a 
certain time period. Such a calculation can be performed using a 
monte carlo type random algorithm in a computer program. 

a 

[0030] Similarly, the movement of the particles can be alterna- 
tively calculated by a monte carlo type simulation of discrete 
particle hopping. For this simulation the hopping probability of a 
particle depends on the potential energy of the initial and the 
final particle position of the hopping event. The potential energy 
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of the particle depends on the material phase. 

[0031] According to an alternative embodiment of the invention, 
the temporal and/or local evolution of the phase interface can be 
described by a movement function of the discontinuity between the 
phases instead of setting up a data field the phase if interface is 
considered to be sharp . 

[0032] The evolution of the phase interface and the movement and/or 
distribution of the particles can be calculated in dependence of a 
high temperature treatment of the material, e.g. an annealing 
treatment . 

[0033] The evolution of the phase interface and/or the movement of 
the particles also can be calculated on the basis of a potential, 
i.e. energetical difference between the two phases. Further, it is 
possible to determine the evolution of the phase interface by 
calculating it in dependence of the concentration of the particles. 
The movement and/or distribution of the particles can be calculated 
in dependence of a driving force induced by the difference of 
potential energy in the two phases. 

[0034] In an alternative embodiment the movement of the particles 
at a discontinuous transition between the two phases can be deter- 
mined in dependence of an accumulation of particles at said 
transition. In this case, a particle transport which is 
proportional to the number of accumulated particles and to the 
velocity of the discontinuous transition between the two phases 
takes place. 

[0035] In a specific case, the method can be applied on a 
semiconductor material, such as silicon, with at least two different 
phases where the particles are dopants, e.g. Arsenic, Boron, 
Aluminum, Indium, Phosphorous or other impurities. 
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[0036] It is possible that one of the phases is amorphous and the 
other is crystalline. Then the method can be applied to determine 
dopant redistribution during SPE. This allows one to study accu- 
rately the dopant redistribution during SPE and calculate effects 
such as Indium dose loss (i.e. movement of Indium atoms out of the 
silicon) , Arsenic dose loss, Fluorine dose loss, Arsenic dopant 
activation, Boron dose loss and Boron uphill diffusion. 

[0037] It is also conceivable that one of the phases is liquid and 
the other is solid or crystalline. The method can e.g. be applied to 
determine dopant redistribution during LPE (liquid phase epitaxy) 
regrowth or to calculate dopant redistribution during float zone 
purification processes . 

[0038] The method can also be applied for determining the dopant 
redistribution during LASER annealing processes where the laser is 
powerful enough to melt or amorphize an illuminated region of the 
semiconductor material 

[0039] According to another aspect of the . invention, a computer 
program product to perform such simulations will be provided, which 
is e.g. usable on a common standard personal computer. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0040] The method will be explained below in a detailed description 
referring to the drawings, wherein: 

Figure 1 is a simulated evolution of the phase interface 
between an amorphous and a crystalline phase of a silicon device 
before, during and after SPE regrowth, and 

Figure 2 shows the dopant distribution in said silicon device 
before, during and after SPE regrowth and the corresponding state of 
the phase interface. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
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[0041] As a practical embodiment of the method, the calculation for 
dopant redistribution during SPE regrowth will be presented in 
detail. It has been implemented in the process simulator FLOOPS 
(FLOOPS User Manual, Release 9.0, ISE Integrated System Engineering 
AG, Zurich, Switzerland, 2003) . In this example the high temperature 
annealing of a silicon wafer is considered, parts of which have been 
amorphized by the implantation of indium atoms at a dose of 10 15 
atoms per cm 2 and an implantation energy of 15 keV. After the ion 
implantation, a temperature treatment of up to 2 minutes at 600°C is 
considered, during which SPE regrowth occurs. 

[0042] The following sequence of two subsequent process steps, 
which can be part of a fabrication process of a silicon device, is 
assumed: 

i. Ion implantation of indium atoms at a dose of 10 15 atoms per 
cm 2 and an implantation energy of 15 keV. 

ii. High temperature annealing of the silicon wafer, which 
leads to recrystallization of the amorphous region. A high 
temperature treatment at 600°C is assumed, which is a typical 
temperature for oxide deposition steps in semiconductor device 
fabrication . 

[0043] The method is used to calculate the final dopant 
distribution after these steps. The procedure for determining the 
final dopant distribution is the following: 

a) Calculation of the distribution of the dopant ions after 
the ion implantation. 

b) Calculation of the crystal damage after ion implantation. 
In particular: Calculation of which part of the silicon is 
amorphous, and which part of the silicon is crystalline. 

c) Calculation of a data field x crystallinity ' , which describes 
the crystal order of the silicon, after the ion implantation. 

In amorphized regions, the value of this data field is 0, in 
perfectly crystalline regions, the value of the data field is 1, in 
crystalline regions with high crystal damage the value is between 0 
and 1 which represents the phase transition region. 
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[0044] These first three steps serve to calculate the data fields 
'dopant concentration' and x crystallinity ' after the ion 
implantation. The next two steps serve to calculate the time 
evolution of these data fields during the SPE regrowth: 

d) Calculation of the time evolution of the data field 
crystallinity during SPE regrowth. The recrystalli zation process is 
reflected in the model by an increase of the value of the 
crystallinity field from 0 (amorphous) to 1 (crystalline) . 

e) Calculation of dopant movement during high temperature 
treatment, including the SPE regrowth. This calculation takes into 
account that a gradient of the crystallinity field can act as a 
driving force for dopant movement. In addition, the usual driving 
forces for dopant movement, as e.g. described in 'P.M.Fahey et al . , 
Point Defects and Dopant Diffusion in Silicon, Review of Modern 
Physics, Vol. 61, No. 2, pp. 289-384, 1989', are also taken into 
account. These include dopant diffusion due to gradients of dopants 
and/or silicon point defects, and dopant drift due to the electric 
field. 

[0045] All five steps serve to calculate the distribution of dopant 
atoms at a certain time, especially the final distribution in the 
steady state at the end of the thermal treatment. During the high 
temperature treatment, the time evolution of the crystallinity and 
of the dopant concentration data fields is calculated by solving a 
set of partial differential equations which describe the time 
evolution of these fields. The partial differential equations are 
solved for a discrete mesh of points in space and time, using 
standard methods of process and device simulation. 

[0046] Below, the above five steps will be explained in more 
detail, with a focus on the steps c-e. 

[0047] The dopant distribution after ion implantation can be 
calculated by major commercial process simulation tools used in the 
simulation of silicon device fabrication, such as DIOS or FLOOPS 
(DIOS User Manual, FLOOPS User Manual, Release 9.0, ISE Integrated 
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Systems Engineering AG, Zurich, Switzerland, 2003) , or it can be 
calculated by monte carlo simulation i.e. tracing a large number of 
ions through the silicon taking into account random collisions with 
silicon atoms and energy loss, or by using analytical formulas, 
which relate the range and distribution of ions in silicon to the 
implantation parameters such as dose, energy, tilt angle. Both 
methods are well-established standard methods. As a result of this 
calculation, one obtains a data field, which describes the as- 
implanted distribution of implanted atoms in the silicon device. 

[0048] The crystal damage is related to the number of silicon atoms 
kicked out from the ideal crystal lattice position to some other 
place in the silicon. It can be calculated by tracing a high number 
of ions on their way through the silicon and count the number of 
collisions which result in a knock-out of silicon atoms, or - 
alternatively - by using analytical formulas, which relate the 
amount and the spatial distribution of the crystal damage to the 
implantation conditions . 

[0049] As for the calculation of the dopant distribution after 
implantation, the calculation of the crystal damage is also a 
standard technology offered by major commercial process simulation 
programs, such as DIOS or FLOOPS (DIOS User Manual, FLOOPS User 
Manual, Release 9.0, ISE Integrated Systems Engineering AG, Zurich, 
Switzerland, 2003) . Therefore the details of this procedure are not 
described here. The crystal damage after ion implantation is stored 
in a data field ^crystal damage' . 

[0050] Regions, where the crystal damage (i.e. the concentration of 
displaced silicon atoms) exceeds a certain threshold, e.g. if more 
than 20% of the silicon atoms are displaced, are considered as 
amorphized. 

[0051] A data field, later on referred as crystallinity , is 
introduced to describe which parts of the device are crystalline and 
which are amorphous. Later, the time evolution of this field is used 
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for calculating the movement, i.e. the local evolution, of the 
interface between the crystalline and the amorphous phase. In this 
example it is assumed, that between perfectly crystalline and 
completely amorphous regions, an extended transition region may 
exist, in which the field changes steadily from 1 to 0 . 

[0052] After ion implantation, the data field crystallization is 
initialized by the following prescription. 

[0053] First, two threshold values are defined for the crystal 
damage from ion implantation. The upper threshold corresponds to a 
displacement of (e.g.) 20% of the silicon atoms. The lower threshold 
corresponds to a displacement of (e.g.) 2% of the silicon atoms. 

[0054] Then, regions, where the calculated crystal damage after ion 
implantation exceeds the upper threshold, are considered amorphous. 
In such regions, the initial value of the crystallization data field 
is set to 0. Regions, where the calculated damage is smaller than 
the lower threshold are considered as crystalline. Here, the initial 
value of the crystallization is set to 1 . In regions, where the 
calculated crystal damage is between the two threshold values, the 
initial value of the crystallinity field changes steadily from 0 to, 
according to the formula: 

crystallization = ( ^crystal damage' - x lower threshold' ) / 
( ^upper threshold' - ^lower threshold' ) 

This means that the crystallization is calculated as the quotient 
between the difference of the crystal damage and the lower threshold 
on the one hand and the difference of the two threshold values on 
the other hand. These regions represent the smooth transition 
between the crystalline and the amorphous part, and typically extend 
over 5-20nm, depending on implantation conditions. 

[0055] According to this embodiment, in the process simulation of a 
high temperature treatment (including diffusion, oxidation, 
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annealing steps), the time evolution of the crystallinity is 
described by an equation. This can e.g. be a partial differential 
equation. As shown in the following, it is possible to formulate the 
equation for the crystallinity field in such a way, that the SPE 
regrowth is well described. 

[0056] Figure 1 shows the calculation of following typical 
situation. For simplicity reasons only one dimension is presented. 
At the beginning of a high temperature treatment a part of the 
semiconductor is crystalline (initial value of the field equals 1), 
and another part of the material is amorphous (initial value of the 
field equals 0) . For the calculation of figure 1 a high tempera- 
ture treatment at 740°C was assumed, which is typical for 
nitride deposition in silicon device fabrication. Starting at 
the initial interface between crystalline and amorphous 
regions, the amorphous region will recrystallize . In the 
model, recrystallization means that the value of the crystal- 
linity is raised from initially 0 to 1 . This recrystallization 
takes place predominantly at the transition region between 
crystalline and amorphous part of the device. Figure 1 shows 
that during SPE regrowth, the transition region is shifted 
towards the amorphous side of the amorphous / crystalline 
interface, and that at the end of the regrowth, the whole 
semiconductor is crystalline. 

[0057] In this embodiment, the following equation for the 
data field crystallinity, which has all desired properties and 
can be solved numerically by a process simulator, is used: 

= div{a grad(C)) + p c\\-C) (C-Y) . (1) 

at 

Here, C represents the crystallinity and is a function of 

space coordinates and time, dC/dt is the time derivative of C, 



15 



div and grad are spatial differential operators, and ot,p and y 
are parameters which may depend on the temperature. In the 
simulation for figure 1, at 740C, the following values were 
used: a=10 _11 cm 2 /s, P=150 s" 1 , y=0.25. These values are chosen 
such that the velocity v SPE of recrystallization (i.e. the 
velocity of the moving interface) equals the value for undoped 
silicon presented in figure 24 of G.L.Olson and J. A. Roth, 
'Kinetics of Solid Phase Crystallization in Amorphous 
Silicon', Mat. Sci . Rep. Vol 3, Nr. 1, June 1988, pp. 1-78. 
The width of the transition region between crystalline and 
amorphous phase is about 15nm. 

[0058] The speed of SPE regrowth depends on the temperature 
and may also depend on the local values of other data fields, 
which are used for the process simulation, such as for example 
the mechanical stress or the concentration of dopant atoms. 
This is taken into account by the choice of the values for a, 
p, and y. a, P, and y can be functions of the temperature and 
of the local dopant concentration. However, in the simulations 
for figure 1 and 2, which are for annealing at constant 
temperature, constant values for a, p, and y were used, 
ignoring a possible influence of the dopant concentration on 
SPE regrowth velocity. 

[0059] Dopant atoms, which are situated in or near the 
transition region between crystalline and amorphous silicon 
are subject to an additional driving force. For dopant atoms, 
which energetically prefer to be in the amorphous part of the 
semiconductor, the driving force will be directed towards the 
amorphous part of the crystalline - amorphous interface. 

[0060] For the origin of such a driving force, one can assume 
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an analogous origin as in the case of the electrical potential 
acting as a driving force for charged particles. The gradient 
of the electrical potential (in other words: the electric 
field) drives positively charged particles towards regions of 
low electrical potential and negatively charged particles 
towards regions of high electrical potential. In the case of 
the crystallinity as the driving force of diffusion, a dopant 
species may have a lower potential energy in amorphous 
regions, where the value of crystallinity is 0, and a higher 
potential in the crystalline region, where the value of the 
crystallinity is 1. As for the electric field, the difference 
in the potential energy causes a driving force, which is 
proportional to the gradient of the potential energy, and 
directed in opposite direction of such a gradient. 

[0061] A possible equation for the driving force F new on 
dopant diffusion due to different values of crystallinity in 
different regions of the semiconductor is: 



F new =-sgrad{C) (2) 

with C being the crystallinity and s being a parameter, which 
depends on the difference of potential energy of an atom in 
the crystalline part and in the amorphous part of the 
semiconductor. Such a new driving force can be applied in the 
process simulation of SPE regrowth. In this example, e=2eV for 
indium atoms was used. 

[0062] The driving force causes a flux of dopant atoms, 
therefore a new term is added to the equation which describes 
the time evolution of the dopant concentration during high 
temperature treatments. The dopant flux if of some dopant A, 
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caused by the driving force F new , can be denoted as : 

<f> = -grad(C)S£A conc (3) 

where 5 is a parameter, which describes the mobility of dopant 
atoms A. As shown in figure 2, A conc is a function of location 
and time. For the calculation of figure 2 (annealing at 600°C 
after indium implantation), the value 8=6. 6- 10" 14 cm 2 /eVs for 
indium atoms was used, and the additional indium flux caused 
by the gradient of the crystallinity was: 



= ~grad{C) - 1 .32 - KT 13 ■ 



[0063] The flux (|) gives rise to an additional term in the 
equation, which describes the time evolution of the 
concentration of the dopant A during high temperature 
treatment : 

dA r 



dt 



= -div{<f>) + other terms = div(grad (C) 5 • e • A conc ) + other terms ( 4 ) 



[0064] In equation (4), 'other terms' represents other 
contributions to the time evolution of A conc , which can be 
independent of the crystallinity field C. E.g., 'other terms' 
can include the conventional diffusion of dopant atoms in 
silicon at high temperatures. div((|>) is the divergence of the 
flux and equals the change of dopant concentration due to the 
flux at each point in the semiconductor. 

[0065] In the example simulation for figure 2, indium 
diffusion was modeled with the Fick' s diffusion law, and it 
was assumed that the indium diffusivity depends on the 
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temperature and on the crystallinity . In this case, 'other 
terms' can be expressed as div ( D(T,C) grad(A conc . ) ), where 
D(T) is the dopant diffusivity, depending on the temperature 
and crystallinity. In this case, equation (4) becomes: 



— f*. = div(grad(C) Se- ) + div(D(T, C) ■ grad(In cnnc . )) ( 5 ) 

at 

where In conc is the indium concentration 

[0066] Figure 2 shows the calculated dopant profile after 
implantation and after SPE regrowth at 600°C, using the 
described procedure. For this calculation the following two 
differential equations were solved with the process simulator 
FLOOPS. 

[0067] Crystallinity equation: 



^ = adiv(grad(C)) + /3C-{\-C)(C-y) . (1) 
at 



With following values at 600°C: a=3-10' 14 cm 2 /s, p=4.8 s~ 1 ,y=0.3 
Equation for indium concentration: 



dln ~~±- = div(grad(C) Se- In COIK _ ) + div(D(T, C) ■ grad{In conc )) 



dt 



With the following parameters at 600°C: 
5 • s = 1 . 32 • 10" 13 cm 2 /s 

D(T,C) = 4.655 • 1CT 18 • ( 1.0 + 3 • 10 7 • ( 1-C ) • C ) cm 2 /s 
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[0068] The diffusivity in crystalline silicon (C=l), 4.655 • 
10" 18 cm 2 /s, is multiplied by a factor, which is 1 in perfectly 
crystalline regions and in perfectly amorphous regions, and 
much larger than 1 in the transition region, where 0<C<1 . 
In this example, a factor of 3.10 7 was used. That means that i 
the simulation model, an increased indium diffusivity is 
assumed in the transition region between amorphous and 
crystalline silicon where the dopants are rearranged from the 
amorphous to the crystalline lattice structure. 

[0069] Figure 2 shows how the dopants have been moved during 
the annealing steps towards the surface of the medium. 

[0070] Equations (1) and (5), or - more general - equations 
(1) and (4), can be integrated into any process simulator, 
which uses continuous fields for the description of dopant 
fields. The equations can be solved with standard numerical 
methods. As the initial distribution of the dopants and the 
initial value of the crystallinity field are known after steps 
1-3 of the procedure, the partial differential equations (1) 
and (4) allow to calculate the final distribution of dopants 
at the end of a treatment, such as a high temperature 
annealing step. 
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